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1. Introduction 

Wilson loops are very fundamental and important quantities in gauge theories. In this 
paper, our main focus will be on their relevance to the description of the infrared (IR) 
behavior of scattering amplitudes. We will consider the general massive case, from which 
the massless one can be obtained as a limit. 

The appearance of Wilson loops in this problem is easy to understand. The infrared 
divergences in scattering amplitudes originate from soft regions of loop integration, for 
which one can employ the eikonal approximation. In this way, one finds that infrared 
divergences of a scattering process are given by a correlation function of Wilson lines, where 
the lines in position space point along the momenta of the scattered particles. However, in 
taking the eikonal limit, additional ultraviolet (UV) divergences are introduced. They are 
equivalent, up to a sign, to the original IR divergences. This allows one to regard the former 



as the UV anomalous dimension of Wilson line operators, whose renormalization properties 
are well understood [1,2]. Note that the latter depends on the color representation of the 
external particles and is in general a matrix in color space. It is known analytically to two 
loops [3] . The analysis of the general structure of this soft anomalous dimension matrix is 
of great interest, with recent studies focussing on the massless case [4,5]. 

In the planar limit, the matrix factorizes into Wilson lines consisting of two segments. 
The cusp anomalous dimension associated to two Wilson lines is known in QCD to two 
loops [6], and in AA = 4 supersymmetric Yang-Mills (SYM) to three loops [7]. In this paper, 
we extend the calculation in planar A/" = 4 SYM to four loops, and, in addition, compute 
the non- planar four-loop value in a special scaling limit. 

The aim of this paper is to develop methods for the computation of such Wilson line cor- 
relators, planar and non-planar, and to deepen the understanding of the functions involved. 
This is closely related to ideas being discussed for understanding the loop corrections to 
scattering amplitudes. The functions that are typically encountered can be described by 
certain classes of iterated integrals. A key problem is to identify which specific class of 
functions is required to describe a given scattering process. It was found that integrals for 
scattering amplitudes or Wilson loops can be put into a 'd-log' form [8-^10], where one can 
pull out an overall normalization factor, and the remaining integrand is a differential form. 
Moreover, such a representation suggests the existence of simple differential equations for 
the integrals. The latter also help to make the transcendent ality properties of the integrals 
manifest. Recently, evidence was presented that integrals having such simple properties 
are not limited to supersymmetric theories, but can be present in generic D = 4 — 2e 
dimensional integrals [11]. The Wilson line integrals considered in this paper can be con- 
sidered as a special, simplifying limit of the more general scattering amplitude integrals. 
We will derive 'd-log' representations for a wide class of Wilson line integrals, relevant to 
the physical problems discussed above, and show how to compute them using differential 
equations. 

The A/" = 4 supersymmetric Yang-Mills (SYM) theory is a good testing ground for 
exploring such Wilson loops. For the specific Wilson loops studied in this paper, results 
can be obtained from various methods such as supersymmetric localization techniques or, 
in the planar case, integrability [12,13]. The AdS/CFT conjecture also allows to compute 
Wilson loops at strong coupling. 

In A/" = 4 SYM, it is natural to define the locally supersymmetric Wilson loop operator 

[14,15] 
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where n is a vector on S^ . It parametrizes the coupling of the Wilson loop to the six scalars 
$. We consider as the integration contour a cusp formed by two segments along directions 
(momenta) p^ and q^, and allow the two segments to couple to the scalars through n and 
n' , see Fig. 1. Then, the vacuum expectation value {W) of the Wilson loop will depend on 



p'^ 




Figure 1: A Wilson line that makes a turn by an angle 4> in Euclidean space. The two segments go 
along p^ and g^, respectively. The vectors n and n' are internal vectors that determine the coupling 
to the six scalars $, see eq. (1.1). 
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as well as on the 't Hooft coupling A = g'^N, and the number of colors N. 

If Auv and Air are short and long distance cutoffs, respectively, then the divergent 
part of the vacuum expectation value of the Wilson loop takes the form [1,2] 



(W) ~ exp 



log Tcusp 

Air 
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(1.3) 



This defines the cusp anomalous dimension rcusp(</'i^i A, A^). 

Note that the dependence of Fcusp on 9 is simple. It can only occur through Wick 
contractions of scalars, and because of SO (6) invariance it appears only through n ■ n' = 
cos 9. Therefore, at L loops, Fcusp is a polynomial in cos^, of maximal degree L. Having 
made this observation, we find it convenient to introduce the variable 



e 
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where the denominator was chosen for future convenience. When the geometric angle <j) 
and internal angle 9 satisfy <j) = i^, which corresponds to ^ — >• 0, the anomalous dimension 
vanishes. Thus we expect the following structure in perturbation theory, 



L = l ^ ^ T = \ 



'■\<t^,l/N^), 



(1.5) 



where A = g^N is the 't Hooft coupling, and g the Yang-Mills coupling. The sum over ^^ 
starts from r = 1, since ^ = corresponds to a supersymmetric configuration, for which 



Fcusp vanishes. 



Note that Fcusp has non-planar corrections starting from four loops. We will discuss 
the full structure of the color dependence to four loops in sections 2 and 4. 



The r = 1 term is known to all loop orders, including the non-planar corrections [16]. 
In this paper, we will compute the full planar result at four loops, for 9 = 0, as well as the 
non-planar contribution to T^'^'^\ which is the leading term in the scaling limit ^ — ?• oo. 
This is done by analytically continuing 6 and keeping i?i> as a free parameter. This scaling 
limit was introduced in ref. [16] and it was shown that it allows to describe the planar 
ladder diagrams in a simple way. This was further developed in [17, 18]. 

This paper is organized as follows. In section 2 we review the color dependence of Fcusp 
to four loops. In section 3 we explain the kinematics and give an overview of the functions 
that will appear in Fcusp- hi section 4 we discuss the structure of a class of Wilson line 
integrals and propose a systematic way of evaluating them. We then apply this to the 
non-planar correction to Fcusp in the scaling limit. In section 5 we give details on how 
we compute the planar Fcusp from a massive scattering amplitude, and give the result for 
the four-loop cusp anomalous dimension in the planar limit. In section 6 we compare our 
results to those at strong coupling. We conclude in section 7. There are two appendices. 
In appendix A we collect the contributing diagrams from the scattering amplitude and 
discuss their structure, while appendix B contains the analytic continuation of Fcusp to 
values beyond threshold. 

2. Color structure to four loops 

Here we discuss the color dependence of Fcusp to four loops. It is best understood using 
results from non-Abelian eikonal exponentiation [19,20]. 

Let us briefly recall our conventions, following [21]. We consider SU(N) with the Lie 
commutator 

^Ta^j.b^^-jabcj.c (2.1) 

and the quadratic Casimir operators 

[T^T%j = CfS.j , r'^f^'' = Ca6''' , (2.2) 

of the fundamental and adjoint representation of the Lie algebra. The fundamental gener- 
ators are normalized by Tr(T'*T*) = TpS""^. Ny^ is the number of generators of the group, 
i.e. the number of gluons. The higher group invariants appearing at four loops can be 
expressed in terms of the following fully symmetrical tensors, 

^abcd r-p Trjiarpb rpC rpd i rparpb rpdrpC i rparpCrpbrpd 

li r? V ix rC tx Ix n, ti, r{, £x ' rC ix n, ti, 

+ nTf^T^T'^ + nTiT'^Tl, + nr^Tf^T'^] . (2.3) 

Here R can be either F or A for the fundamental and adjoint representation, respectively, 
with [T^]fec = [^"]6c and [T^]fec = —if°'^'^- For SU(N) generators in the standard normal- 
ization we have 

2 ' ' 27V ' ^ ' Na 48 ^ ^ 



From [19,20] it follows that Abelian-like terms containing powers of Cp cancel in Tcusp; 
thanks to the logarithm in its definition, see eq. (1.3). Moreover, an analysis of the possible 
color diagrams reveals that the result for Fcusp at one, two, and three loops is proportional 
to CftCfCatCfCJ^, respectively. At four loops, we can have two structures, CpC"^, and 
the quartic Casimir d'^'^'^d'^'^'^/N. 

In summary, we have, to four loops 

log{W) = g^Cpwi + g^CFCAW2 + g^CpClw^ + 5' 



jabcdjabcd 



WAb 



, (2.5) 



where we have chosen the normalization (W) = 1 + 0{g'^) . Here the webs Wi correspond 
to linear combinations of Feynman diagrams. The explicit expressions are easily obtained 
by the method of [19,20]. 

One advantage of this formulation is that one can directly compute that logarithm of 
the Wilson loop correlator, and that each web only has an overall divergence^. The latter 
is easy to remove, so that in practice one can define Fcusp in terms of finite integrals. 

Using eq. (2.4), we make the dependence on A^ of eq. (2.5) manifest. 
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(2.6) 



We see that in the large A^ limit, keeping A = g'^N fixed, only the contribution g^N^W4b 
disappears from the R.H.S. of eq. (2.6). In other words, to three loops, it is sufficient 
to know the planar result for the Wilson loop in order to restore the full result in eq. 
(2.5). At four loops, an additional computation of the diagrams contributing to w^b is 
required. In the remainder of this paper, we compute the non-planar contribution Wib in 
a recently-introduced scaling limit, as well as the full planar result to four loops. 

3. Kinematics and integral functions 

As explained in the introduction, we will mainly be interested in the (p dependence of Fcusp. 
Here we discuss a convenient kinematical variable, and different physical regions. We also 
introduce a class of functions that we find appropriate to express the answer in, and discuss 
the branch cut structure of the latter. 



3.1 Kinematical structure 

It is convenient to introduce a new variable x = e**^, which in general is complex. The com- 
putation we are considering is invariant under (/) — t- —(/). This corresponds to an inversion 
symmetry in x. 

There are three different kinematical regions that we would like to discuss. It is useful 
to recall the relationship of Fcusp to IR divergences of scattering processes involving massive 



^We tacitly assume that the intrinsic renormalization of the bare parameters of the Lagrangian has 
aheady been carried out. 



particles, such as e~^{p) — )■ 7*e^(g), which have the same analytical structure. (See e.g. 
refs. [22,23].) With the on-shell conditions p^ = q^ = rn? (in the mostly-minus metric 
H — •••—), this process is naturally described using the variable s/m?^ where s = {p — qf'. 
It is related to x via 



y/j - A7V?/S - 1 

X = ^= . (3.1) 

V^l - 4:w?/s + 1 

There we distinguish three kinematical regions, above threshold s > Am?, below threshold, 
< s < 4m^, and finally the space-like region s < 0. They correspond to regions III, I , 
and II, respectively, that we now discuss from the Wilson loop viewpoint. 

Region I: The first region corresponds to real (j), (f) £ [0,it]. This means that the 
absolute value of x is 1. In this case we have a cusp in Euclidean space, and Tcusp is real. 
The two limiting cases are the following: for (p = the contour is a straight line, and Tcusp 
vanishes (for 6 = 0). The first correction ~ 0^ in this small angle limit is known exactly 
in A and N [16,24]. The opposite limit </) — )■ vr is related to the quark-antiquark potential. 
This limit is subtle and requires a resummation of certain diagrams, see [7,17,25-27]. One 
may also view the Wilson loop as the eikonal approximation to a form factor of massive 
quarks. In that case, this region corresponds to the region below the two-particle threshold. 

Region II: We can analytically continue (f> to Minkowskian angles. In that case, x is real 
and positive. Because of the inversion symmetry x — t- 1/x, it is sufficient to take x £ [0, 1]. 
The second endpoint, x = 1, again corresponds to the case of a straight line discussed 
above. Near the endpoint x — t- 0, on the other hand, the cusp anomalous dimension 
diverges linearly in log(x), to all orders in the coupling constant [6]. The coefficient of the 
linear divergence is the well-studied light-like cusp anomalous dimension; the latter can 
also be obtained from the anomalous dimension of high spin operators [28-30] . We may 
remark that the Wilson loop approach considered here is a very efficient way of computing 
this quantity. 

Region III: Finally, we have the region above the threshold of creating two massive 
particles. This region corresponds to x being real and negative. As before, it suffices to 
take x £ [—1,0], because of the inversion symmetry in x. Tcusp has a branch cut along 
the negative real axis, and the iO prescription in the propagators implies that x has a 
small imaginary part. For the mostly minus Minkowski-space metric that we are using, the 
position-space propagator connecting two segments of the Wilson loop is proportional to 
(s^ + t'^ + st{x + l/x) — iO)^^^*^, where the line parameters s and t are positive, and hence 
for X G [—1,0] we should add a small positive imaginary part to x. In this region, Tcusp 
has an imaginary part. 

3.2 Harmonic polylogarithms and symbols 

What are the functions needed to describe Tcusp? Results at lower loop orders and for 
related scattering processes suggest that the class of functions we are seeking are the har- 
monic polylogarithms (HPL) [31]. They are generalizations of ordinary polylogarithms, 
and appear naturally e.g. within the differential equation technique to evaluate loop inte- 
grals, see e.g. [22,23]. They are also natural from the point of view of the singularity and 



branch cut structure described in the previous paragraph, with x = 0, ±1 being special 
points. They are defined iteratively by 

Ha,,a,,...,aSx)= I faAt)Ha,,...,a„it) dt , {oi, 02, . . . a„} / {0, 0, . . . 0} , (3.2) 

Jo 
where the integration kernels are defined as 

h{x) = {l-xr\ fo{x) = x-\ /_i(x) = (l + x)-i. (3.3) 

The degree (or weight) 1 functions needed to start the recursion are defined as 

Hi{x) = -log(l - x) , Ho{x) = log(x) , H^iix) = log(l + x) . (3.4) 

There is a special case when all indices are zero, Hq _ _ ^ q(x) = ^ log"(2;). The subscript 




of H is called the weight vector. A common abbreviation is to replace occurrences of m 
zeros to the left of ±1 by ±{m + 1). For example, Hofi,ifi,-i{x) = H^-2{x). 

HPLs have simple properties under certain argument transformations, and one can 
use their algebraic properties to make their asymptotic behavior manifest. We refer the 
interested reader to ref. [31]. A very useful computer algebraic implementation has been 
given in refs. [32,33]. For fast numerical evaluation, especially at complex arguments of 
the HPLs, we found the C++ implementation in GiNaC [34] very helpful. 

As we will describe in the section 5, we compute the planar Fcusp from a massive 
scattering amplitude, where at each loop order a certain number of individual integrals 
appears. It turns out that each of these integrals can be expressed as a linear combination 
of HPLs of argument x where in general all possible weight vectors appear at a given degree. 
In the total result, however, we find the simplification that the result can be written in 
a compact form when using HPLs of argument 1 — x^, and weight vectors with indices 
0, 1 only! The latter property is also present in the four-loop non-planar correction in the 
scaling limit, and becomes manifest from the formulas in sections 4 and 5.3. 

In the context of the iterated integrals and differential equations studied in section 
4, the notion of the symbol of a transcendental function [35-37] is very useful. It can 
be derived recursively for any function fwixi, ■ ■ ■ ,Xn) of weight w whose total differential 
assumes the form 

dfw = '^fi,w-idlogRi, (3.5) 

i 

where the fi^w-i are of weight w — 1 and the Ri are algebraic functions. The symbol S{fw) 
is then defined recursively via 

S{U) = ^S{fi,^^i) <S) R^, (3.6) 

i 

which involves a tensor product over the group of algebraic functions. We emphasize that 
eqns. (3.5) and (3.6) make the close connection between the 'd-log'-representations (to be 





Figure 2: Propagator- type integral discussed in the main text. 

discussed in the next section) and the symbol of a function manifest. We also note that 
symbols of the HPLs discussed above are built from the alphabet {x, 1 it x}. As a specific 
example, we have 

S{Hn{x)) = cS(Li„(a;)) = -(1 - x) (g) x (g) • • • (^ x . (3.7) 

^^ ^ ' 

n— 1 terms 

4. Iterated Wilson line integrals, and nonplanar result in scaling limit 

Here we discuss a general method for computing Wilson line integrals in position space. 
We then apply it to the computation of the non-planar cusp anomalous dimension in the 
scaling limit. 

4.1 'd-log' forms for Wilson line integrals 

In this section we elaborate on 'd-log' forms for integrals, which were introduced in the 
context of scattering amplitudes in [8-10]. As an instructive example, let us discuss the 
diagram shown in Fig 2. The corresponding integral over the line parameters s and t can 
be written as 

/■ dS Adt X f ,, , N , n / ^ / . - N 

/ -^ 7^ 77 TT^ = -. / dlog(s + tx) A(ilog(t + sx), (4.1) 

A s^ + t^ + st{x + 1/x) 1 - x2 7^ ^^ 

where on the RHS we have dropped differentials involving dx because they do not contribute 
to the integral, and where the integration region A is s G [a, S] and t G [b, T]. 

What is gained from writing the integral in this way? We see that a natural normaliza- 
tion factor, x/(l — x^), has been pulled out of the integral. Together with trivial prefactors 
originating from the Feynman rules, this constitutes the normalization of the diagram. The 
remaining integral will give a (generalized) polylogarithmic function, which, in the present 
example, has degree 2. It depends on the variables a, b, S, T and x, 

/(a, b, S,T,x) = / d log(s + tx) A d log{t + sx) . (4.2) 

Ja 



Integrals of this type satisfy simple differential equations, as we explain below. Let us first 
focus on one of the two integration variables, say s, and rewrite the integral in a more 
convenient form thanks to the identity 



s ~\~ (y 
d log(s + a) ^ d\og{s + 13) = d log A d log(a - /3) . (4.3) 

s ~\~ p 

A simple generalization of this identity holds for n-forms. Then, we perform one integration 
at a time, in this case starting with the one over s. The main point is that one will always 
have an integral of the form 



Gia,f3i):= d\og{y + a)F{y, (3i) , (4.4) 

JAy 

where y is the integration variable, and a and /3j are parameters, and F is some function. 
Then the algorithm outlined in appendix A of ref. [38] can be used to determine the 
differential of G. It can be expressed in terms of quantities appearing in the differential of 
F. In our example, a short calculation gives 

^ ^^ I. o rr \ Ji u^ {b + ax){S + bx) {T + ax){a + bx) 

'^(°'^'^'^'"^=^^°g^^°g (a + 6x)(b + 5x) +'^°g°^°g (5 + ax)(a + rx) 

{b + Sx){S + Tx) {T + Sx)ia + Tx) 

+ "^°g^^°g (s + fex)(r + 5.) +'^°g^^°g (g + r.)(r + ax) 

+ dlogxlog (^ + "")(^ + ^"^^^ + ^"^(" + ^") (4 5) 

This equation determines / up to an integration constant. The latter can be determined, 
for example, from the boundary condition that / vanishes at a = 5. Note that it is trivial 
to read off the symbol of / from eq. (4.5). 

In the present example, one can also directly integrate eq. (4.5). The answer obtained 
can be written in terms of dilogarithms, 

f{a,b,S,T,x) =Li2 [-g^j ~^^'^[~Sx) ~ ^^^ ( ^) ^ ^^^ \~~ax 

-L„(-|.)+L,.(-|i)+Li,(-^.)-Li,(-^i) . (4.6) 

This agrees with ref. [39]. 

In summary, we see that one can always compute the symbol of integrals that are of the 
type (4.1), and for generalizations with more propagators, and poly logarithmic functions 
inserted into the integrand. In particular, any ladder integral appearing in Tcusp can be 
computed in this way. Preliminary results suggest that the generalization to graphs with 
interaction vertices is possible. For example, in ref. [18], for two classes of diagrams the 
internal integration associated to the interaction vertex could be computed analytically, 
with the remaining integral of the form (4.4). 

We used this method to compute the non-planar ladder integrals appearing in Tcusp 
to four loops. In the next subsection, we discuss which integrals are required, and in the 
following subsection the results are reported. 






(a) 



(b) 



(c) 



Figure 3: All diagrams contributing to the quartic Casimir terms at four loops, in the scaling 
limit. 

4.2 Nonplanar contribution to scaling limit 

Here we compute the mtegrals contributing to w^h of eq. (2.5) in the scaling limit. Thanks 
to the scaling limit, we only need to keep ladder diagrams with four rungs between the two 
Wilson line segments. Moreover, only diagrams containing the color factor d'^'^'^d'^'^'^ are 
required.^ 

It is easy to see that the only ladder type diagrams containing the quartic Casimir are 
the ones shown in Fig. 3. Denoting their color factors by Ci and the integral functions by 
Xj, with i = a,b,c we have 



log(l^). 



Ca Xa + Cblb + "^ Cc 1c + 



(4.7) 



where the dots represent other diagrams with color factors involving only Cp and Ca, 
and the 2 is a combinatorial factor, due to the fact that Fig. 3(c) also appears with the 
two Wilson lines interchanged. The color factors of these diagrams contain a trace over 
eight generators, e.g. Tr(rT''r^r'^TT^r'=T'^) in Fig. 3(a), and similarly for the other two 
diagrams. Using the Lie commutator one can show that up to terms proportional to powers 
of Cf and Ca, Tr(T"T^T=T'^r'^r''r^r'^) is given by Tr(T"r^r^T^)Tr(r|T^r|rj), which 
in turn is related to d'^'^'^d'^'^'^, see table 11 of [21]. Exphcitly, we have 

1 



lahcd labcd 



arj-idrr-iC rjib ^ 



Tr(r"r''r^r'*)Tr(riT^rir| 



12 



NCfC\ , 



rr\ /rT-iarpbrpCrj-idrj-iarpbrpCrT-id'S 



dfcd^fcd ^ ^^^ 



Cp — 3CpCA H — T P' A 



11 

T 



19 

24 



cl 



(4.8) 
(4.9) 



Moreover, using the Lie commutator (2.1) one sees that the color factors of the three 
diagrams in Fig. 3 are related. One finds (we normalize all color factors by Tr[li7'] = N), 



Cc=Ca + -CpCa{Cf - CA/'i){Cp - Ca) , 
Cc=Ca- -CpCa{Cf - Ca/'2) . 



(4.10) 
(4.11) 



^Recall that the color dependence of a general four-loop diagram can be expressed in terms of Cf, Ca, 
and d'^'"^d'^^'"^/N. The terms with powers of Cf higher than one cancel in log(W), as per eq. (2.5). 



10 



This means that the three diagrams contribute equally to the color factor dp (F^ /N . 
Taking this into account, we find that the term proportional to d'^'^'^d'^'^'^/N is given by 
(in the scaling limit), 

W4b^la+lb + 21c. (4.12) 

Let us discuss the definition of the integrals. They are line integrals of the type considered 
in section 4.1. Here a comment on the regularization of the Wilson loop operator is due. 
Naively, it has both infrared as well as ultraviolet divergences. We are interested in the 
ultraviolet divergences. Tcusp is defined as the coefficient of the ultraviolet divergence. 
Since log{W) only has an overall divergence, it is easy to see how different regularization 
procedures are related. The position space calculations above can be formulated e.g. in 
cut-off regularization for both IR and UV divergences. Another possibility is to treat the 
integrals as in heavy-quark effective theory (HQET), with dimensional regularization. 

In both cases, one can make the logarithmic divergence transparent by changing vari- 
ables. Let us denote the line parameters on the two lines by Si and tj, with i = 1, ... 4, 
respectively. After rescaling all variables Sj = psi,ti = pii, with X]i=i(^« + ^«) — 1> ^^^ P 
integral contains the divergence. When working with cutoffs, this integral takes the form 

— =logT ■ (4.13) 

Auv P ^UV 

On the other hand, in HQET with dimensional regularization, one obtains 

^e- = i^ + 0(.°) . (4.14) 

In both cases, the coefficient of the p integral is the contribution to the cusp anomalous 
dimension that we wish to compute, and it is given by convergent integrals. 

For concreteness, let us choose the cutoff version of the calculation. In that case, taking 
into account the discussion above and writing the integrals in d-log form as in section 4.1, 
we have 

li = log -^ ( ]^ _ ^2 ) ^»' i = a,b,c. (4.15) 

Next, we can use the algorithm of section 4.1 to derive iterative differential equations for 
Xj. From these equations, we can immediately determine the symbol of these functions as 
a corollary. We find that they are given by symbols composed from the alphabet x, 1 — x^. 
This implies that they can be expressed in terms of a subset of the HPLs discussed in 
section 3, namely those with indices drawn from 0,1, if we choose x^ or 1 — x'^ as argument 
of the HPLs. 

In order to determine the full functions from the differential equations, we have to 
complement them by boundary conditions. The kinematical point x = 1, or equivalently 
(/) = 0, is a good boundary condition, where the Feynman integrals are expected to be 
regular. However, the prefactor (x/(l — x^))^ in (4.15) diverges in this limit, and hence the 
functions /j must have corresponding zeros. We find it likely that a careful investigation 
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of the iterative differential equations would reveal that this boundary condition fixes all 
undetermined constants. We found that simply using the condition of regularity of eq. 
(4.15) at a; = 1 determined most coefficients, and we computed the remaining ones by 
considering asymptotic limits x — t- 0, which we evaluated using standard Mellin-Barnes 
techniques. For more details on such methods, see section 5. In this way, we found 



Xa = - 6tT^Hi^1,1,2 + 48i7i,i,i,4 - 8^2^1,1,2,1 + 64i/i,i,2,3 + 64ffi,i,3,2 

— 67r2i7i,2,i,i + 48i7i,2,l,3 + 48i/i,2,2,2 - 107r2i/i,i,i,i,i + 80i7i,i,i,i,3 
+ 80ifi,i,i,2,2 + 24ifi,i,i,3,i + 647/1,1,2,1,2 + 32i:fi,i,2,2,l + 32/^1, 1,3,1,1 

+ 48/^1,2,1,1,2 + 24i/i,2, 1,2,1 + 24f/'i,2,2,i,i + 62ffi, 1,1,1,1,2 + 40ffi, 1,1, 1,2,1 

+ 22//i,i,i,2,i,i + 8-^1, 1,2, 1,1,1 + 6//i,2,i,i,i,i + -ffi, 1,1, 1,1, 1,1 , (4.16) 

1 f 

ib = - 47r2i?i,i,i,2 - y 7r2//i,i,2,i + 16i7i,i,2,3 + 32i/i,i,3,2 - ^TT^Hi,2,l,l 

20 a 

+ 16//l,2,l,3 + 16/^^1,2,2,2 - Y^^i/1,1,1,1,1 + 16^1,1,1,1,3 + 24/^-1,1,1,2,2 
+ 24//i,i,2,l,2 + 8i/i, 1,2,2,1 + 16i/l,l,3,l,l + 16-f/'i,2,l,l,2 + 8/^1,2,1,2,1 

+ 8i/i,2,2,i,i + 40/^1,1,1,1,1,2 + 24ifi, 1,1, 1,2,1 — 8i/i, 1,2,1, 1,1 — 8i/i,2,i,i,i,i 

+ 4ifi,i,i,i,i,i,i + 48C3i^i,i,i,i , (4.17) 

X, = + 47r2i/i,i,i,2 - 12i/i,i,i,4 + y Vr2i7i,i,2,i - 28i/i,i,2,3 - 40//l,l,3,2 

+ 47r2//i,2,i,i - 24i/i,2,l,3 - 24i/i,2,2,2 + YVr2i/i,i,i,i,i - 32/71,1,1,1,3 

— 38//i,i,i,2,2 — 6i/l, 1,1,3,1 — 34//i,i,2,l,2 — 147/1, 1,2, 2,1 — 20ifi,i,3,i,i 

— 24//i,2,l,l,2 — 12i/l,2, 1,2,1 — 12//l,2,2,l,l — 38//i, 1,1, 1,1,2 — 22//i, 1,1, 1,2,1 

— 4i/i,i,i,2,i,i + 2i/i, 1,2,1, 1,1 + 2i/i,2,i,i,i,i + 2i/i,i,i,i,i,i,i — 24^37/1,1,1,1 . (4.18) 

Here we use the abbreviation H^ = Hw{l — x^). Recalling eqs. (4.15) and (4.12), this 
determines w^i, in the scaling limit. 



5. Planar calculation from massive scattering amplitude 

It was shown in ref. [40] that the Regge limit s/m'^ ^ 1 of the planar Coulomb branch 
amplitude M{s/m'^,t/m'^) is governed by the cusp anomalous dimension. This connection 
was used in [7] to compute the three- loop value of Tcusp- 

The advantage of this approach is that an expression for the integrand of M{s/m?,t/m?) 
is already known. It is in the form of a small number of scalar integrals, each of which 
results from many Feynman diagrams. This simple integrand was obtained by using gen- 
eralized unitarity, in conjunction with (extended) dual conformal symmetry [41,42]. It has 
been pointed out that the limit relating the amplitude and Fcusp also works at the level of 
the integrand [18] . This implies that one can obtain an efficient integral representation for 
Fcusp in this way. 
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Here we wish to extend the work of [7] to four loops and determine the planar part 
of Fcusp from the four-loop scattering amplitude. As a starting point, convenient Mellin- 
Barnes representations for the eight contributing integrals are available from [42]. The 
strategy of our calculation is the following: First, we use generalized cuts to determine 
the power of ^ to which each scattering amplitude integral contributes. Details of this 
procedure can be found in appendix A. Next, based on experience from lower loop orders 
and the structure observed for the ladder diagrams [18], we make an ansatz subject to 
certain assumptions, which we will specify below. This reduces the calculation to the 
determination of a certain number of undetermined coefficients. In order to determine the 
latter, we analyze both the Mellin-Barnes representations and the ansatz in various limits, 
such as j; — )• and j; — )• 1. In this way, we obtain (more than) enough algebraic equations 
to determine all coefficients. Moreover, this provides consistency checks for the ansatz. 

5.1 Assumptions and classification of HPLs 

An analysis of the cuts of the integrals contributing to the scattering amplitude suggests 
that the four-loop result for 6 = 0, where we have ^ = (1 — x)/(l-|-x), has the structure 

rcusplAV(8.^)^ = E (r^)''r^'^''^(^) + o{i/N') , (5.1) 

where the T^^''^'{x) are certain transcendental functions. What can we assume about their 
structure? Looking at the results up to three loops we may make a number of observations. 

• all results for Tcusp can be written in terms of harmonic polylogarithms 

• the degree of transcendentality of the functions involved is uniformly (2L — 1), where 
L is the loop order 

• the subset of HPLs with indices 0, 1 only and argument \ — x^ is sufficient to describe 
the answer. In terms of the symbol, this means that only letters x, 1 — x^ are required. 

In the case of ladder integrals, the first two items can be proved, and the third item is true 
at least up to six loops [18]. As we showed in section 4, it is also true in the non-planar 
case. We find it reasonable to assume that these properties hold true for the full result at 
four loops. 

If this assumption is correct, the calculation is reduced to the determination of the 
precise linear combination of the allowed functions. Our starting point will be all functions 
Hw{l — x'^) of weight seven, where the weight vector w is build from entries 0, 1. We also 
allow transcendental constants such as Q, QCji possibly multiplying lower degree functions 
to construct a term of total degree seven. 

We can restrict and classify these functions according to their symmetry properties. 
In fact, Fcusp has to be symmetric under the inversion x — t- 1/x. This follows from the 
definition cose/) = l/2(a; -|- l/x). As a result, r(^'^)(a;) has to be odd/even under this 
symmetry for r odd/even. In this way, we find 51 even and 50 odd functions under this 
symmetry. 
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Another simple condition we can impose is that T^'^'^'{x) should have at least r zeros 
as X — >• 1. The reason is that it is multiplied by the factor ^'", which for 9 = 7r/2 has a 
degree r pole at rr = 1. But x = 1 corresponds to the straight line case, and this should be 
finite for each integral contributing to Fcusp- We will also verify this behavior by expanding 
integrals near x — t- 1. 

In order to further restrict the number of functions, we make two additional observa- 
tions about the results at lower loops L < 3, eqs. (5.4)-(5.9). Inspecting them one sees that 
in these cases the number of O's in the weight vector of the functions T^ '"^'(x) is always 
smaller than r. We will assume this to hold true also at four loops. Moreover, T^ '"^'(x) 
vanishes for L>latx = — 1.^ This is required in order to obtain the correct leading order 
behavior at rr — )■ — 1, which corresponds to the quark antiquark potential limit. This limit 
will be discussed in more detail at the end of this section. 

Imposing all these conditions our ansatz becomes, in summary, 

• 12 functions for r(^'^)(x), of degree 7, indices 0, 1, even under x — )• 1/x, at most one 
entry in weight vector, and vanishing as (1 — x)^ as x — )• 1 and as (1 + x)"'^ as x — )■ —1. 

• 21 functions for r'^'^'(x), of degree 7, indices 0, 1, odd under x — t- 1/x, at most two 
entries in weight vector, and vanishing as (1 — x)^ as x — t- 1 and as (1 + x)^ as 
X — > -1. 

The terms F^^'^' and F^^''*-' were already computed in refs. [16] and [18], respectively. 

5.2 Asymptotic limit of Mellin-Barnes integrals 

Let us now explain how to determine the coefficients of the ansatz. By means of the 
Mathematica packages MB.m [43] and MBasymptotics.m [44] we perform the asymptotic 
expansions of the Mellin-Barnes representations, first about the point x = 0. The expansion 
parameter appears in the form x^ log'^(x), and each of these terms is accompanied by 
one or several transcendental constants. For p > these constants are in general not of 
homogeneous weight, but the highest transcendentality is always 7 — q."^ We determine 
these constants analytically for p = 0, . . . , 6 and q = 3, . . . ,7. For p = we also include 
q = 2. It is interesting to note that at most two-dimensional Mellin-Barnes integrals are 
required for this calculation at g > 2, and three-dimensional ones at q = 2. After we 
computed all relevant constants in this way, we perform the series expansion of our ansatz 
to the respective orders in x and log(x), and solve the resulting algebraic equations for the 
unknown coefficients appearing in our ansatz. As an illustrating example, take 



22 27r2 
¥"^"3" 



X log°(x) = ( -;^ ai + -^ a2 ) + TT^ ( -- ai + - 02 



r/ 8 28 \ 




-—-ai + --a2 


+ TT" 


[V 27 9 J 





X log^(x) , (5.2) 



where the LHS stems from the solution of the Mellin-Barnes integrals at a particular power 
in X and log(x), and the RHS stands for the expansion of the ansatz to the same order. 

^In this case one has to rewrite r^'^''^'{x) in terms of HPLs of argument x. We will discuss the analytic 
continuation of Fcusp in the next section and in appendix B. 
^Except for q — 6, where it is zero. 
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Assuming 1 and vr^ to be linearly independent we obtain two algebraic equations, yielding 
ai = —3 and 02 = 1/2. 

Since some coefficients in our ansatz appear only at powers q = 1 or q = 0, the above 
procedure does yield most, but not all coefficients. In order to determine the remaining 
ones we expand the Mellin-Barnes representations about the point x = 1. In this limit the 
expansion is purely of the form (x — 1)'^, without any logarithms. We include all terms 
with s < 4 and determine the coefficients in the same way as above. However, this time 
we have to solve Mellin-Barnes integrals that are up to seven-dimensional. 

We emphasize that the number of equations this procedure yields is much larger than 
the number of undetermined coefficients in our ansatz, such that our system of equations 
is largely overconstrained. As a rule of thumb we have about twice as many equations 
compared to the number of coefficients. This property is of utmost importance since 
otherwise potential inconsistencies in our ansatz could not be revealed. 

We also do numerical checks, but only after the application of MBasymptotics.m, i.e. 
we check numerically all analogues of the x-independent part of the LHS of eq. (5.2). 
Performing numerical checks on the unexpanded expressions is not well suited here since 
by construction the integral and the ansatz agree to high powers in x and {x — 1). Hence 
the ansatz obtained in this way will agree very well numerically with the integral we are 
computing, even if the ansatz was incomplete. 

Last but not least we have the algebraic cross-check that the final answer does only 
diverge linearly in log(x) as x — t- 0, see eq. (5.14). This cross-check is non-trivial since it 
connects different powers of ^, each of which diverges with the seventh power of log(x). 

5.3 Planar result to four loops 

We are now in the position to present the results for the cusp anomalous dimension up to 
four loops in the planar limit. We have 

4 / X \ L L 



1j — 1 T — 1 



^)(x) + 0(A^l/iV2), (5.3) 



where 

r(i^i)=^/7i, (5.4) 

r(^^^) = - ^^1,1,1 - \-'H, , (5.5) 

r(^^^) =\h,,, + \h,,,,, (5.6) 

at one and two loops [6,45-47], 

r^'^^) =\^'H,,,,i + ^^1,1.1,1,1 + ^^1 , (5.7) 

r^^'') = - ^C3^1,l - \^^Hi^2 - \t^^H2,1 - -yHi^l,l - ^1,1,1.2 - 7^1,2,1.1 

lb 6 4 4 ' 

— -f^2, 1,1,1 — ^-f^i,i,i,i,i , (5.8) 
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113 

r*- ' =Hii^s + -H'1,2,2 + -^1,1,1,2 + --ffi, 1,2,1 + --ffi,2,i,i + T-ffi,i,i,i,i (5.9) 

at three loops [7], and 

r^^''^ = - l^'^i,i,i - ^'^1,1,1,1,1 - 1^1.1,1,1,1,1,1 - ^^'i^i , (5.10) 

5 2 45 

4'i 2 1 'i 13 

r(^'2) =^C5^1,l + gVr^Cs^l,! + 5(3^1,1,1,1 + Y^^'^1,2 + Y^^'^2,1 + ^^'^1,1,1 

+ TT -f^l,l,l,2 + 71" -f^l, 1,2,1 + tTT -f/'l,2,l,l + ^VT i:f2,l,l,l + T^VT i:/'l,l,l,l,l 

TO O ^ 

+ 5i:/'i_i,i,i,i,2 + --ffl, 1,1, 1,2,1 + ;^-f^l,l,l,2,l,l + 3.^1,1,2,1,1,1 + -^.f^l,2,l,l,l,l 

+ -r-f^2, 1,1, 1,1,1 + -r7r-f^i,i,i,i,i,i,i ) (5.11) 

4 lb 

r(4;3) ^ _ 3^3//, ,2 _ 4C3ii'l,2,l - 3C3^2,1,1 - 5(3^^1,1,1,1 " J^ ^'^1,1,1 " ^^'^1,1,3 

^ O r> Q Q I Q "^ O 

— o^ -^1.2,2 - TT .^1,3,1 - TT .f^2,l,2 " VT ii'2,2,1 " -TT .f/'l,l,l,2 " -TT i:^l, 1,2,1 

On Q ^y Q / 

— 77 TT .f/^1,2,1,1 — VT .^2,1, 1,1 — TTT^ .f^l,l,l,l,l — 5ffi 1 1 1 3 — -i:/^l,l,l,2,2 
24 2 

9 

— 3ii'i,i,2,l,2 — 2i:/'l,l,2,2,l — 3/^1,1,3,1,1 — 5/^1,2,1,1,2 — 4/^1,2,1,2,1 — :^//l,2,2,l,l 

43 

— 3//i,3, 1,1,1 — 5/^2, 1,1, 1,2 — 6//2, 1,1,2,1 — 5/^2,1,2,1,1 — 3//2,2, 1,1,1 ^//l, 1,1, 1,1,2 

;^//l, 1,1, 1,2,1 — 8//i,i,i,2,i,i — 7//1, 1,2, 1,1,1 ^//i, 2,1, 1,1,1 ;^//2, 1,1, 1,1,1 

— ^//i,i,i,i,i,i,i , (5.12) 

lo 

r^ ' =3//i, 1,1,4 + 4//1, 1,2,3 + 4//1, 1,3,2 + 3//i,2,l,3 + 3//l,2,2,2 + 5//l,l,l,l,3 + 5//i,i,i,2,2 

3 3 

+ 7^ //l, 1,1, 3,1 + 4//l,l,2,l,2 + 2//1, 1,2,2,1 + 2//1, 1,3,1,1 + 3//i,2,l,l,2 + 7j//l,2,l,2,l 

3 23 

+ 7^//l,2,2,l,l + ^//l, 1,1, 1,1,2 + 5//1, 1,1, 1,2,1 + 4//1, 1,1,2, 1,1 + 3//1, 1,2, 1,1,1 

9 23 

+ t//i,2,i,i,i,i + ^//i, 1,1, 1,1, 1,1 (5.13) 

at four loops. As mentioned above, the terms r^^'^^ and r*-^'^-* were already computed in 
refs. [16] and [18], respectively. The remaining terms F'^'^' and F^^''^^ are new. We derived 
them analytically, subject to the assumptions discussed in the previous section. 

In the above equations H^ := //^(l — x^), and x = e**^. The perturbative results given 
in this section and in section 4 can be straightforwardly evaluated numerically in the region 
II, i.e. < a; < 1. Other regions can be reached by analytical continuation, respecting the 
branch cut properties discussed in section 3. We collect the relevant formulas in appendix B. 

A curious feature of the result up to four loops, already remarked upon in [18], is that 
once the result is written in terms of HPLs with argument \ — x^ as above, all HPLs in 
(^_l)('"+i)p{i;'") come with non-negative coefficients. 

In Figures 4 and 5 we plot the cusp anomalous dimension in Regions II and I, respec- 
tively. From the plots one can see the properties discussed below and in section 3. 
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Figure 4: Plot of (— l)^"'"^rcusp as a function of (/3 = — log(a:) with x G [0, 1], i.e. Region II. From 
bottom to top the plot shows L = 1,2,3,4. The small and large (p behavior is known to all loop 
orders: For small (p the first term is quadratic, with the coefhcient given by the Bremsstrahlung 
function. At large (p, Fcusp grows linearly, with the coefficient determined by the light-like cusp 
anomalous dimension. 



We now consider various limits of Fcusp- First, we can use the above results to analyt- 
ically compute the light-like cusp anomalous dimension. It is obtained by taking the limit 
X — ;■ 0, where Fcusp diverges logarithmically, 

1 



lim F, 



cusp 



log(x) F°° + ^0 + 0{x) 



(5.14) 



We find 
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(5.15) 



This agrees with previous numerical results at four loops [41,42,48,49], and with the spin 
chain prediction from ref. [50]. The behaviour (5.14) can also be seen from Fig. 4, where 
the curves grow linearly for large values of (/9 = — log(a;). For Qq we find 
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Go is related to the collinear anomalous dimension for mass-regulated scattering ampli- 
tudes [40]. Unlike F°°, this quantity depends on the regularization scheme and takes a 
different value in dimensional regularization [51]. 
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Figure 5: The functions (— l)^rcusp in the interval (f> € [0,7r], i.e. Region I. From bottom to top 
the plot shows L = 1,2, 3, 4. 



This behaviour can be seen from Figs. 4 and 5, where all curves start quadratically from 
the respective origin. The expansion in (5.17) is in perfect agreement with the four-loop 
expansion of the exact result in [16,24]. 



The third limit to consider is x 
and (5 — )• 0. We find 
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with Ls = ln(25/e). Figure 5 shows the divergences as (/> — t- vr. 
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The limit x — >■ — 1 is related to the quark antiquark potential. This limit is subtle. 
Due to ultrasoft effects, a resummation is required. This is done by matching fixed order 
calculations against an effective field theory calculation. In the context of AA = 4 SYM, 
this was discussed in ref. [26], and more recently in [27]. 

We close this section with two remarks. First, note that for 6 = 0, we have ^ = 
(1 — x)/(l + x). It is natural to expect that the full 6 dependence can be obtained by 
replacing ((1 — x)/(l + x)Y in eq. (5.3) by ^'', see section 1. 

Our second remark concerns the regularization scheme dependence. The above method 
assumed a supersymmetric regularization scheme, and therefore we expect our result for 
Tcusp to be valid in that scheme. The transition to other schemes, such as MS, is discussed 
in ref. [52], and has been explicitly worked out there to two loops. 

6. Comparison to strong coupling via AdS/CFT 

Can we compare the fixed order perturbative results of section 5.3 to the results available 
at strong coupling via the AdS/CFT correspondence? The authors of ref. [53] proposed 
such a procedure in the case of the light-like cusp anomalous dimension. They combined 
perturbative data with the string theory insight that the strong coupling expansion takes 
the form 

Fcusp = cVX + ... , (6.1) 

where c is negative, and we work in the planar limit. In order to incorporate this behavior 
they proposed the following ansatz /(A) for Fcusp 

2n 

A- = ^a[/(A)]\ (6.2) 

r=n 

where n is connected to the loop order L via n = L — 1. The constants Cr can be fixed using 
perturbative information. Of course, one can also use strong coupling data, as in [41] , in 
order to gain insights on weak coupling. Here we will use the perturbative two-, three-, and 
four-loop results for Fcusp in order to determine the coefficient of vA at strong coupling. 

Let us give more details about this ansatz in the simplest case, i.e. n = 1. Here the 
extrapolation is based on the two-loop perturbative information, i.e. 

Fcusp = ^^lA + V2X^ + 0{X^) . (6.3) 

Then one can determine the coefficients in eq. (6.2) to be Ci = l/vi,C2 = —V2/V1. The 
latter equation then implies that the ansatz for the interpolation function is 



^«-£ 



-1 + W1-4A- 

Vl 



(6.4) 



which at strong coupling gives "vAv ~:^ -|- . . .. This procedure can be generalized to 
higher loops, where eq. (6.2) implies that one has to solve equations of higher degree, 
which can be done numerically. 
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We then compare this extrapolation to results obtained via the AdS/CFT correspon- 
dence. At strong coupling, i.e. A » 1, Wilson loops are described by minimal surfaces [14]. 
We have the expansion (6.1), and focus on the case 6 = 0. The result for c for region I, i.e. 
real angles, is available through implicit equations involving Elliptic integrals from [54]. In 
the case of region II, we can use the formulas of [55]. The first subleading coefficient in the 
expansion (6.1) is also known [56]. Let us discuss the results of the comparison for the two 
regions in turn. 

Extrapolation for region I: We used the ansatz of eq. (6.2) in order to extrapolate the 
strong-coupling coefficient c in eq. (6.1) from the knowledge of the four-loop data. We 
found that the extrapolations based on our four loop results and the vA behavior give 
a leading order strong coupling answer that agrees to within 2 per cent for the range of 
(j) G [0.1,2.5]. However, for > 2.5 the relative error grows significantly. This is not 
surprising since there one approaches the quark-antiquark limit (^ — )• vr. 

Extrapolation for region II: Here we find very good agreement between the extrap- 
olation based on the four-loop perturbative data, and the strong coupling answer. It is 
interesting to note that the relative error to the strong coupling value goes down from 
approximately 25%, 3% and under 1.6%, when using two-loop, three-loop and four-loop 
data as input, respectively. It is also remarkable that this relative error stays small for all 
data points analyzed in the interval x S [0, 1], despite the fact that the leading coefficient 
(in front of vA) at strong coupling varies by several orders of magnitude. 

Let us comment on the radius of convergence of the expansions. It is known for 
the X —7- 1 and x — t- limits, respectively. The former is described by the Bremsstrahlung 
function [16], whose perturbative series has a radius of convergence of Ac ~ 14.7. The latter 
is governed by the light-like cusp anomalous dimension, where the radius of convergence [50] 
is Ac = vr^. 

7. Conclusion and Outlook 

We computed the velocity-dependent cusp anomalous dimension in maximally supersym- 
metric Yang-Mills theory to four- loop order. The result can be expressed in terms of 
harmonic polylogarithms of degree seven, with argument 1 — x^ and non-negative indices 
only. We determine the non-planar correction at four loops in the scaling limit, which 
involves quartic Casimir invariants as color factors. The method of 'd-log'-representations 
for iterated Wilson line integrals turns out to be extremely powerful for this purpose. It 
allows one to compute the symbol of such functions. If the symbols correspond to a known 
class of functions, HPLs in our case, one can integrate back using boundary conditions. 

Moreover, we determine the full planar four-loop result from massive scattering ampli- 
tudes, where we use asymptotic expansions of Mellin-Barnes integrals to analytically pin 
down the coefficients of a well-motivated ansatz. Our analytical result gives the correct 
values of the four-loop light-like cusp anomalous dimension that was previously calculated 
only numerically [41,42,48,49]. 

We also compare our perturbative result to strong coupling, and find that our extrap- 
olation agrees to better than two per cent with the corresponding string theory result, over 
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a wide range of parameters. 

Taken together, the only pieces missing to obtain the full - planar and non-planar - 
result of the velocity-dependent cusp anomalous dimension to four loops are the non-planar 
terms proportional to ^^ and ^^. The light-like limit of the non-planar cusp anomalous 
dimension is also envisaged in [57] by means of the on-shell form factor. However, we 
emphasize that the present approach allows to obtain the full j;-dependence, and not just 
the light-like limit. 

The results we have derived here shed light on the structure of the planar four-particle 
amplitude on the Coulomb branch of 7\^ = 4 super Yang-Mills. The latter is an infrared 
finite function M{s/'w?' ^t/rn^). Kinematically it is very similar to light-by-light scattering 
via massive particles. It is an interesting open question what class of two- variable func- 
tions describe such processes beyond the one- loop order. The integrals we have computed 
determine the asymptotic limit of this amplitude as s/m? ^ 1. 

There are several generalizations to the 'd-log'-approach discussed in section 4. The 
first generalization concerns the Wilson loop contour. While we have focused on a contour 
formed by two segments in this paper, it is clear that the technique applies equally to 
contours formed by n segments meeting in a point. This is relevant for the description of 
infrared divergences of massive scattering amplitudes at the non-planar level, see e.g. [3,58]. 
We also wish to emphasize that massless results can be obtained as a corollary. 

Another obvious generalization of the applicability of this technique has to do with 
the regularization. On physical grounds, at least in principle, one can always choose com- 
binations of diagrams that only have a superficial UV divergence. For such quantities, one 
can easily switch between regulators. Our method is very naturally formulated in a cut-off 
scheme, however it is equally possible to use dimensional regularization. This is straightfor- 
ward for integrals that only have a superficial UV divergence. For other integrals, one first 
has to identify the integration regions that lead to divergences and perform subtractions. 

Finally, whereas we focussed in the present paper on scalar and guon exchanges, pre- 
liminary results suggest that the generalization to graphs with interaction vertices is pos- 
sible [18]. 
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A. Four-loop integrals and generalized cuts / leading singularities 

The velocity-dependent cusp anomalous dimension TcuspC^) can be obtained from the Regge 
limit of massive amplitudes in A/" = 4 super Yang-Mills [40,42]. At the four-loop level, there 
are eight contributing diagrams, which are depicted in Figure 1 of [42]. The corresponding 
amplitude is given by eq. (2.8) of that reference. In the Regge limit s — >• oo the logarithm 
of the amplitude is given by^ 

logA^ '^ log(-mVs) reusp(-mVt) , (A.l) 

where t is related to x via —rri^/t = x/(l — a;)^. 

We now study the systematics of the Regge limit at the four-loop level [42] by consid- 
ering the integrals contributing to the four-loop amplitude. We expect them to have the 
general structure 

/ = /ox/, (A.2) 

where /q is an algebraic normalization factor, and / is a function having degree of tran- 
scendentality eight. In the literature, such functions are sometimes referred to as pure 
functions. 

Generalized cuts or leading singularities are useful in order to test whether (A.2) holds, 
and to determine the normalization factor Iq. 

For example, for the massive box integral at one loop, normalized by st, we have 



/o~ l/Vl-4m2/s-4m2/t- (A.3) 

Notice that in the Regge limit, this factor becomes proportional to ^. 

Likewise, we computed the maximal cuts of all integrals up to four loops. The result 
is consistent with eq. (A.2), and we find the following behavior of the prefactors as s — >• oo 
(the superscript "r" means that the integral is rotated, i.e. s ■v^ t) 



ha ~ C, 


iia ~ e 


hh^ i , 


II, - e 


he ~ e, 


-'4c ~ ? 


hd ~ e, 




he ~ C, 


rr t2 
-'46 ~? 


hf ~ i, 


rr i-2 


hd2 ~ i, 


rr e2 

hd2 ~ ? 


hf2 ~ c^ 


-'4/2 ~ ? 



(A.4) 

Notice that due to exponentiation, the maximal power of Regge logarithms that a given 
integral has is bounded by the power of S,- Comparing to Appendix A of [42], we find that 
this is in agreement with the above ^ dependence. 

We see that we can classify the contribution of the integrals to Fcusp according to 
which power of ^ that they are normalized by. This is a very useful feature, as it allows to 
compute the contributions to different powers of (, independently. 
^Note that refs. [40,42] use difTerent metric conventions. 
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B. Analytic continuation of Fcusp 

In order to analytically continue Fcusp to regions I and III it is sufficient to apply the 
argument transformation 1 — x^ — )• x^ to the HPLs in section 5.3 and subsequently extract 
the logarithms explicitly. This gives 

r(i;i)=-log(rE), (B.l) 

r(2;i)= + llog3(:,) + !^log(x), (B.2) 

r(2;2) = _ 1 log3(^) _ ^ log(x) - log(rE) H2{x^) + H^{x^) - C3 (B.3) 

3 6 

at one and two loops, 

r(3;i) = _ 1 iog5(^) _ ^ iog\x) - ^ log(x) , (B.4) 

b 6 b 

r(3;2) = _ 6i75(x2) + hog\x)H^{x^) - ^log3(x)F2(x2) +'^-^\og\x)H^{x'') 

vr^ 9 

- \og^{x)H^[x^) - — \og[x)H2{x^) + - \og{x)Hi{x^) + Cs log^lx) 

+ — log5(x) + -^ log3(x) + -^ log(x) + 6C5 , (B.5) 

3U 9 3d 

2 

r(3;3) = _ 2(:^H2{x') + ^H^{x^) + ^H^{x^) + 2H2,3{x^) + 3H-i,2ix^) + 3H^,iix^) 

b 

- hog^{x)H2{x^) + log\x)Hs{x^) -'^log{x)H2{x^) -2log{x)Hi{x^) 

1 27r^ 

- 21og(x)if2,2(2;2) - 21og(x)if3,i(a;^) - Cs^og'^i^) - -log^lx) - --log3(x) 

5 9 

-|)logM-?fe-^C3 (B.6) 

at three loops, and 

r(^^^) = + ^ log^(x) + ^ log5(x) + ^ log3(x) + ^ log(x) , (B.7) 

r(4;2) = - ^log7(x) - ^H,{x')log^{x) + ^-H2{x')logHx) - ^vr2log5(x) 

- ^7r^Hi{x')log\x) - 2H,{x^)log\x) - Calog^Cx) + ^^^^sC^^) log3(x) 

+ ^^4(x2)log3(x) - i^vr4log3(x) - ^7r4ifi(x2)log2(x) - ^-7r'H3ix^)log\x) 
+ yi?5(x')log'(x) - yC5log'(x) - ^^2^3log'(x) + ^7r^i72(:E2)log(x) 
+ ^^2i/4(x^)log(x) - ^F6(x2)log(x) - ^log(x) - In^Hsix') 

- In'H.ix') + 2.f>H,{x^) - 35C7 + ^ + ^ > (B-S) 
239 38 1 4 

r(^^=^) = + — log^(x) + ^^l(x') l0g6(x) - -H2{X^) l0g5(x) + -TT^ log5(x) 
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y 

- ^7r2i72(x2) log3(x) + mAix") log3(x) - ^i?l,3(x2) log3(x) - 2F2,2(x2) log3(x) 

4 47 17 

-4F3,i(x2)log3(x) + -i7i(x2)C3log3(x) + -— 7r4log3(x) + -7r4i7i(x2) log2(x) 
6 1U8 45 



2^ 
3' 



+ ^7r2i73(x2) log2(x) - 10F5(x2) log2(x) + 27r2i7i,2(x2) log2(x) + 9i74(x2)C3 



+ 27r^i72,i(x^)log^(x) - 2i/2,3(a;')log'(a;) - 2/^3,2 (x^) log' (x) + 8(5 log' (x) 

+ 2F2(x2)C3log2(x) + ^7r2C3log2(x) + l^T:^H2{x^)\og{x) - -Ti^ H^{x^)\og{x) 
+ 20i?6(x') log(x) - 27r2i:fi,3(x2) log(x) - 2i7i,5(x2) log(x) - ^'K^H2,2{x') log(x) 

+ 7F2,4(x2) log(x) - ^7r2F3,i(x2) log(x) + 10i/3,3(x') log(x) + 16i:f4,2(a:') log(x) 
+ 22/75,1(2:') log(x) + 2/7i(x2)C5log(x) + 2C|log(x) + 2^2^i(x2)C3log(x) 

- 4i/3(x')C3log(x) + 7r6log(x) - -rzTT^H^{x^) - 32Hrix^) + -tt^ H2,3ix'') 

945 45 3 

- 10H2,5ix^) + 7r2i73,2(x2) - 13/73,4(x') + tt"^ H^^iix^) - ISH^^six^) - 28H54x^) 

- 40if6,i(x2) + I6C7 + lOH2{x^)C5 + vCs - ^vr2i72(x2)C3 + ^ , (B.9) 

2 6 45 

^(4^4) ^ _ 46 j^^7 _ 3^^ 2) j^g5(^) _ 23^2 i^g5(^) ^ ^^(^2) i^g4(^) _ ^^ i^g4(^) 

315 5 9U 

- ^7r2i72(x2)log3(x) - ^i74(x2)log3(x) - 2i72,2(:r2) log3(x) - ^//3,i(a;')log'(x) 

19 2 

- — ^4iog3(x) + -7r2i/3(x2) log2(x) + iH.ix^) log2(x) + 3i72,3(a;') log'(x) 

+ 6//3,2(x2)log2(a:) + 9i74,i(x2)log'(x)-^C5log'(x)-3/72(x2)C3log'(x) 

- ^7r2C3log2(x) - ^7r4i72(x2)log(x) - ^7r2F4(x2)log(x) - AHq{x^) log{x) 

- Tr^H2,2{x^) log(x) - 6H2,4{x^) log(x) - -Tr^H3,i{x^) log(x) - 10i73,3(x2) log(x) 

- 16i?4,2(x2) log(x) - 22i75,i(x2) log(x) - 6H2,2Mx^) log(a^) - 6i:^2,3,i(a:') log(x) 

- 8H3,i,2{x^) log(x) - 8i:f3,2,i(a;2) log(x) - 6/^4,1,1(2:=') log(x) - Cl log(x) 



177r'' 



1 



1 



+ 2Hs{x^)C3 log(x) - — - log(x) + -n^H^ix^) + -n^H^ix^) + 8Hr{x^) 
252U 3U 2 

+ ^7:^H2,3{x^) + 9i72,5(x') + vr2i73,2(x2) + 14i73,4(x2) + ^n^H^^^ix^) 

+ 19i/4,3(x2) + 25i:f5,2(2;2) + 30/76,1(2:') + QH2,2Mx^) + 9i^2,3,2(x2) + 9H2,4Ax^) 
+ 8H3,i,3{x^) + 14i:f3,2,2(x2) + 14if3,3,l(2:') + 15i:f4,i,2(x2) + 15H4,2,l{x^) 



+ 12/7,i,i(x^)-^J 



H2{x^)C5 - ^ - Jvr2i72(x2)C3 - 4H4x^)C3 



TT 



QH2,2ix')C3 - 8H3,lixnC3 - ^CS 



30 



(B.IO) 
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at four loops. In region III, i.e. x G [—1, 0], the logarithms are the only source of imaginary 
parts. Together with the iO-prescription from section 3 the imaginary part can therefore 
be extracted explicitly in this region. 
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